arXiv: 1504.0351 lv2 [hep-ex] 14 Sep 2015 


EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN) 




CERN-PH-EP/2013-037 

2015/09/15 


CMS-SMP-13-013 


Measurement of the Z boson differential cross section in 
transverse momentum and rapidity in proton-proton 

collisions at 8 TeV 


The CMS Collaborations 


Abstract 

We present a measurement of the Z boson differential cross section in rapidity and 
transverse momentum using a data sample of pp collision events at a centre-of-mass 
energy y/s = 8 TeV, corresponding to an integrated luminosity of 19.7 fb The Z 
boson is identified via its decay to a pair of muons. The measurement provides a 
precision test of quantum chromodynamics over a large region of phase space. In 
addition, due to the small experimental uncertainties in the measurement the data 
has the potential to constrain the gluon parton distribution function in the kinematic 
regime important for Higgs boson production via gluon fusion. The results agree 
with the next-to-next-to-leading-order predictions computed with the FEWZ program. 
The results are also compared to the commonly used leading-order MadGraph and 
next-to-leading-order POWHEG generators. 
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1 Introduction 

The production of lepton pairs in proton-proton collisions is dominated by the Drell-Yan (DY) 
process i.e. the production of an intermediate 7*/Z boson by the incoming partons. Mea¬ 
surements of the cross sections as a function of the mass of the intermediate boson (hereafter 
referred to as the 'Z boson'), rapidity, and transverse momentum provide a very sensitive test 
of quantum chromodynamics (QCD). Precise measurement of the differential cross section also 
allows comparisons to calculations employing different parton distribution functions (PDF) 
and underlying theoretical models. Finally, the understanding of DY lepton pair production 
is important in the study of several physics processes, such as diboson and tt production, as 
well as in searches for new resonances decaying to dileptons in models of physics beyond the 
standard model. Differential measurements of Z boson production at the LHC have already 
been performed UJ-jB]. 

In this Letter we present the first measurement of the DY cross section at a centre-of-mass en¬ 
ergy of 8 TeV for dimuon pairs in the vicinity of the Z boson peak, doubly differential in the 
transverse momentum qj and in the rapidity y of the Z boson. The analysis uses the data sam¬ 
ple of pp collisions collected with the CMS detector at the LHC in 2012, corresponding to an 
integrated luminosity of 19.7 fb 1 . We present the absolute fiducial cross section and the fidu¬ 
cial cross section normalised to the inclusive fiducial cross section. The measurement probes 
the production of Z bosons up to high transverse momenta, qj ~ 100 GeV, a kinematic regime 
where the production is dominated by gluon-quark fusion. The precision of this measurement 
leads to experimental uncertainties smaller than or similar to the uncertainties of the gluon 
PDF in the kinematic region that is relevant to the production of the Higgs boson via the gluon 
fusion mechanism. Using the Z boson production process to constrain the gluon PDF |9) in 
the future would be complementary to other processes such as direct photon production llOl 
and top-quark pair production fill that constrain the gluon PDF in this regime. Moreover, sev¬ 
eral of the experimental systematic uncertainties in the DY measurement are uncorrelated with 
these other processes. The latter have more complex topologies and thus have complementary 
and potentially larger systematic uncertainties. 


2 The CMS detector 

A more detailed description of the CMS detector, together with a definition of the coordinate 
system and the relevant kinematic variables, can be found in Ref. [12l . The central feature of 
the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a mag¬ 
netic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead 
tungsten crystal electromagnetic calorimeter (ECAL), and a brass/scintillator hadron calori¬ 
meter (HCAL), each composed of a barrel and two endcap sections. Muons are measured in 
gas-ionisation detectors embedded in the steel flux-return yoke outside the solenoid. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 
Muons are measured in the pseudorapidity range |?/| < 2.4, with detection planes made using 
three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. Matching 
muons to tracks measured in the silicon tracker results in a relative px resolution of 1.3-2.0% 
in the barrel and better than 6% in the endcaps, for muons with 20 < pj < 100 GeV. The px 
resolution in the barrel is better than 10% for muons with px up to 1 TeV [13). The particle-flow 
event reconstruction [141 115] is used in this analysis. It works by reconstructing and identify¬ 
ing each particle with an optimised combination of all subdetector information. The energy 
of photons is directly obtained from the ECAL measurement, corrected for zero-suppression 
effects. The energy of electrons is determined from a combination of the track momentum at 
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the main interaction vertex, the corresponding ECAL cluster energy, and the energy sum of 
all bremsstrahlung photons. The energy of muons is obtained from the corresponding track 
momentum. The energy of charged hadrons is determined from a combination of the track 
momentum and the corresponding ECAL and HCAL energies, corrected for zero-suppression 
effects, and calibrated for the nonlinear response of the calorimeters. Finally, the energy of neu¬ 
tral hadrons is obtained from the corresponding calibrated ECAL and HCAL energies. The first 
level of the CMS trigger system, composed of custom hardware processors, uses information 
from the calorimeters and muon detectors to select the most interesting events in a fixed time 
interval of less than 4 ps. The high-level trigger processor farm further decreases the event rate 
from around 100 kHz to around 400 Hz before data storage. 

3 Simulation 

The signal process is simulated using the leading-order (LO) MadGraph 1.3.30 IT16I generator 
with 0-4 additional jets, interfaced with PYTHIA HTTl v6.4.24 with the Z2* tune Itl8l . The match¬ 
ing between matrix element calculation and parton shower is performed with the /<x-MLM al¬ 
gorithm Ifl9l . Multiple-parton interactions are accounted for via PYTHIA. The LO CTEQ6L1 Il20l 
PDF set is used for the generation. As a cross-check, a second signal sample is simulated using 
the next-to-leading-order (NLO) POWHEG ll2Tl - [24l generator interfaced with PYTHIA. For this 
generation the NLO CT10 |25l PDF set is used. 

The backgrounds are generated with MadGraph (W+jets, tt, rr), POWHEG (single top quark 
1261127 ]), and PYTHIA (dibosons, WW, WZ, ZZ). The inclusive cross sections of DY, W+jets ||28| , 
and tt Il29l processes are normalised to next-to-next-to-leading-order (NNLO) predictions. In 
addition, for the single top quark a higher-order (approximate NNLO Il30l l inclusive cross 
section is used. The generated events are passed through a detector simulation based on 
GEANT4 13TI . The simulated processes are overlaid with minimum bias collisions in order 
to reproduce the distribution of the number of additional proton-proton interactions per bunch 
crossing (pileup) present in data. 

In addition, for comparison with the final result the double differential cross section is com¬ 
puted with fewz 3.1.b2 Ii32l at NNLO. The electroweak corrections are computed at NLO and 
initial-state photon radiation and photon-induced processes are included in the generation. 
The computation is done for each cjj bin separately. The factorisation and renormalisation 
scales are chosen as V M| + q where Mz is the mass of the Z boson and qj is the value of the 
lower edge of the corresponding bin in qj. For the computation the NNLO NNPDF23 PDF set 
with radiative corrections ff33l is used. 

4 Event selection 

An isolated single-muon trigger is used with a threshold of px > 24 GeV and a requirement of 
|;/| < 2.1. The standard CMS baseline offline muon selection Ifl3l is applied. It requires that 
the muon candidate is reconstructed both in the muon detectors and in the inner tracker, with 
X 2 / Hdof < 10 for the track fit. In addition, requirements are placed on the minimum number 
of pixel and tracker layers that are hit and on detailed matching criteria between the trajecto¬ 
ries reconstructed in the inner tracker and the muon system. The distance between the muon 
candidate trajectory and the primary vertex is required to be smaller than 2 mm in the trans¬ 
verse plane and smaller than 5 mm in the longitudinal direction. The vertex with the highest 
sum of pj of associated tracks is selected as the primary vertex. The leading reconstructed 
muon in px is required to be the one selected by the trigger. In order to be within the trigger 
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acceptance, the leading muon is selected with pj > 25GeV and | ?/1 < 2.1. The second muon is 
required to have pj > 10 GeV and )/] < 2.4. The relative isolation is defined to be the scalar 
sum of the transverse momenta of charged hadrons, neutral hadrons, and photons in a cone 
of A R = V (A;/) 2 + (A cp) 2 < 0.4 around the muon direction, divided by pj. After correction 
for pileup, the value of the relative isolation is required to be less than 0.12 (0.5) for the lead¬ 
ing (second) muon. A pair of oppositely charged muons is required to have an invariant mass 
M(pp) between 81 and 101 GeV. In the rare case of ambiguity among several reconstructed 
muons, the muon pair with the invariant mass closest to the Z boson mass is selected. The 
absolute rapidity \y\ of the muon pair must be less than 2 . 

Scale factors are applied to account for known differences between data and simulation. The 
efficiencies for the tracking, the trigger, the muon isolation and identification are determined 
via a "tag-and-probe" method Il34l . The tracking efficiency is measured in bins of //. The trigger 
efficiency is measured in bins of muon pj and rj for positive and negative muons separately. 
The identification efficiency is measured in bins of pj and //. In particular phase space re¬ 
gions, especially for higher qj, the second muon can often point opposite to the Z boson in 
the azimuthal plane. In that direction, the hadronic activity from the recoil of the Z boson is 
enhanced. Thus the second muon is often less isolated than the leading muon and the isolation 
depends on the event kinematics. For that reason, the requirement for the isolation of the sec¬ 
ond muon is looser than the requirement for the leading muon, and the efficiency is measured 
in variables reflecting the second muon direction with respect to the Z boson. Three variables 
for the second muon are chosen to measure this effect on the efficiency in data: the transverse 
momentum of the dimuon system qj, the cosine of the polar angle cos 9* and the azimuthal 
angle <p*. The two angles are measured in the Z boson rest frame, where the z axis is the Z 
boson flight direction. For cos(0' ) = — 1 the leptons are more likely to be close to the hadronic 
recoil. The azimuthal angle is chosen to be zero for the proton closest to the z axis in this frame. 
The isolation efficiency for the leading muon is measured in bins of pj and //. These efficiencies 
are measured in data and simulation, and scale factors are applied to the simulation to account 
for differences with respect to the data. 

The backgrounds are small relative to the signal (at the percent level or smaller) and can be 
divided into two categories: those where the leptons come from Z boson decays and those 
where the leptons stem from other sources. The backgrounds from tt, TT, WW, tW, and W+jets 
are estimated from specific data samples. Backgrounds typically have two prompt leptons, al¬ 
though not necessarily of the same lepton flavor: flavor universality is used for the background 
estimation. The estimation consists of two steps. First, the oppositely charged mixed lepton e /7 
yields are measured in both data and MC. Then the ratio of the yields in data and simulation 
in this data sample (e /7 channel) is used to normalise the simulation in the muon channel. The 
e/7-channel selection uses the same trigger as the final sample, thus the same trigger efficiency 
scale factor is used. In addition, the tracking, the identification, and the isolation efficiency scale 
factors for the leading muon are applied. Electrons are selected if they have pj(e) > 20 GeV 
and 1 77 (e) | < 2.1, which is similar to the fiducial regions of the muon selection. No data-to- 
simulation correction factors are applied to the electron identification since the effect on the 
final results is negligible. In order to enhance the statistical precision, the invariant mass range 
of the e /7 pairs is increased to [60,120] GeV. Within the uncertainties, no significant trend in qj 
and | y\ is observed in the ratio of the e/7 yields in data and simulation, and a constant scale 
factor of 0.987 ± 0.008 is used. The WZ and ZZ backgrounds, which include a true Z —> / 7/7 
decay, are taken from simulation. 
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5 Measured observables and granularity 

The reconstructed and background-corrected double differential distribution in qj and \y\ is 
unfolded to pre-final-state radiation (FSR) lepton kinematics. The unfolding is performed to 
the kinematic region 81 < M(yy) < 101 GeV and within the kinematic selection of the leading 
(second) muon, pj > 25(10) GeV and |?/| < 2.1 (2.4). The unfolding is done using an iterative 
unfolding technique lf35l implemented in the RooUnfold package It36l . The bins in qj are [0,20], 
[20,40], [40,60], [60,80], [80,100], [100,120], [120,140], [140,170], [170,200], [200,1000], In \y\ a 
constant bin width of 0.4 is used and the binning ends at 2. 

MadGraph is used as simulation input to the unfolding. The unfolding is validated by treat¬ 
ing the simulated POWHEG signal sample as data. The unfolded POWHEG distribution is found 
to be compatible with the distribution at the generator level within unfolding uncertainty. 


6 Systematic uncertainties 

The sources of systematic uncertainty are ordered by their average size starting with the largest 
one. The full covariance matrix is computed for both the normalised and the absolute cross 
section. 

• Luminosity uncertainty: 

The uncertainty in the measurement of the integrated luminosity is 2.6% 11371 . 

• Tracking, muon trigger, isolation, and identification efficiency correction factors: 

A potential bias in the measurement of the efficiencies with the tag-and-probe tech¬ 
nique is estimated by varying the most sensitive components: the background in 
simulation is removed and doubled; the signal is parametrised with the sum of two 
Voigtian functions instead of the sum of a Crystal Ball and a Gaussian function; the 
efficiencies are parametrised only in ?/ but with finer bins; and, only tags with a sin¬ 
gle available probe are selected for the measurement instead of all possible pairs. For 
each contribution a 100% correlation is assumed in the covariance matrix. The effect 
of statistical uncertainties in the measured data-to-simulation scale factors is esti¬ 
mated by their variation within the uncertainties in a series of pseudo-experiments. 
Combining the effects extracted from these variations, the systematic uncertainties 
are typically between 1% and 1.6%, depending on the bin, and increase with qj. 

• Pileup uncertainty: 

The cross section of minimum bias events is varied by ±5% and the impact of the 
pileup multiplicity in the simulation on the measurement is used as correlated un¬ 
certainty for all bins. This uncertainty is at maximum around 0.5% and is negligible 
compared to the leading uncertainties. 

• Statistical uncertainties of the simulation: 

The uncertainty due to the limited number of events in simulation is estimated via 
pseudo-experiments by varying the response matrix and the efficiency within the 
statistical uncertainties. 

• FSR: 

The simulation is reweighted to reflect the difference between a soft-collinear ap¬ 
proach and the exact O(a) result, similar to what was done in Ref. Il34ll . It also reflects 
effects from higher-order contributions. The difference between the measurements 
with and without the reweighting is assigned as an uncertainty and is assumed to 
be fully correlated for the covariance matrix. 



5 


• Backgrounds: 

• tt,tW, WW, W+jets, and rr backgrounds: 

A10% uncertainty is assigned to the scale factor derived in the ey method. 

This accounts for the statistical uncertainty of the scale factor and for the 
uncertainties in the lepton efficiencies. For the covariance matrix full cor¬ 
relation is assumed. 

• WZ and ZZ backgrounds: 

The diboson backgrounds that include a Z boson are determined from 
simulation. The cross sections are varied by 50% to estimate the system¬ 
atic uncertainty. While the inclusive cross sections have been measured 
to agree reasonably well Il38l - l40ll . we assign conservatively 50% to account 
for the fact that we use the qj and \y\ shapes from LO calculations. 

• Muon momentum resolution: 

The muon momentum resolution is measured in data and simulation, and corre¬ 
sponding corrections are applied. The covariance accounting for the statistical un¬ 
certainty of the muon momentum correction measurements is calculated via pseudo¬ 
experiments. In addition, an uncertainty is assigned to take into account possible 
correlated offsets. 

• Z boson polarisation: 

The lepton angular distribution of the Drell-Yan process can be described at LO 
through the coefficients, Aq-A^ Hill . However, inaccuracies in the way this is mod¬ 
elled in the simulation can affect the result of the unfolding. The angular coefficients 
A 0 -A 4 are inferred in bins of c]j and \y\ in Il42ll for both data and simulation we use. 

For each parameter Aj the simulation is independently reweighted to correspond to 
the data as measured in l42l . In case the difference in Aj is smaller than the typical 
theoretical uncertainty of 10% [|43l Aj is varied by 10%. The full difference between 
the default polarisation and the changed polarisation is assigned as systematic un¬ 
certainty. Full correlation is assumed between the bins. 

• qj and y shapes: 

The dependence of the results on the qj and y shapes of the simulation is studied by 
repeating the analysis using POWHEG as the signal sample. The results obtained us¬ 
ing MadGraph or POWHEG for the unfolding are compatible with each other within 
the statistical uncertainties. In addition, the MadGraph simulation is weighted in 
fine bins in qj and y to match the background-corrected data. The effect on the 
result using the reweighted simulation for the unfolding is much smaller than the 
uncertainties assigned to the limited statistics of simulation and is neglected. 

The contributions of the uncertainties to the normalised cross section measurement are pre¬ 
sented in Fig. |lj The systematic uncertainty is dominated by the uncertainty in the efficiency 
correction. In the highest bins of qj the measurement is dominated by the statistical uncertain¬ 
ties. The uncertainty contributions to the absolute cross section measurement are presented in 

Kg-0 

7 Results 

The double differential cross section normalised to the inclusive cross section for Z bosons 
decaying to muons is presented in Table [l] A comparison of the measurement with the NNLO 
FEWZ computation is shown in Fig. |3j where the first five plots show the i]j dependence in the 
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Figure 1: Relative uncertainties in percent of the normalised fiducial cross section measure¬ 
ment. Each plot shows the qj dependence in the indicated ranges of \y\. 
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Figure 2: Relative uncertainties in percent of the absolute fiducial cross section measurement. 
The 2.6% uncertainty in the luminosity is not included. Each plot shows the qj dependence in 
the indicated ranges of \y\. 


five bins in \y\ and the last plot shows the qj dependence integrated over \y\. In the bottom 
panels the ratio of the FEWZ prediction to data is shown. The vertical error bars represent the 
statistical uncertainties of data and simulation. The red-hatched bands drawn at the points 
represent the systematic uncertainties of the measurement only. The scale uncertainties are 
indicated by the grey-shaded areas and the PDF uncertainties by the light-hatched bands. The 
scale uncertainties are estimated from the envelope of the following combinations of variations 
of the factorisation u p and the renormalisation //r scales: (2}i^,2y^), (0.5/p,0.5//k), (2 
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(0.5piF,}iR), and (ji[--,0.5 fi[{). The PDF uncertainties are evaluated as the envelope of the 
uncertainties of the NNLO NNPDF23 and the NNLO CT10 If44l PDF sets. The scale uncertainty 
is about 4% for the lowest qj bin. In the second q-y bin it is about 8% and increases up to about 
14% in the highest qj bin. The jump in the size of the scale uncertainty between the first and the 
second bins in qj can be understood as a consequence of reducing the order of the calculation 
to NLO when the Z boson is produced in combination with a jet, which is the dominant process 
for qj > 20 GeV. While the scale uncertainties are smaller at low qj, the shape is not expected 
to match the data well since multiple soft gluon emissions are not modelled. At very high qj 
QED corrections could reach a few percent [45,46]. 

The PDF uncertainties in the region qj > 20 GeV range between +1% and —4%. The uncer¬ 
tainty is asymmetric because the inclusive cross section computed using the NNLO CT10 PDF 
set is about 2.5% larger than the one obtained using the NNLO NNPDF23 PDF set. 

The NNLO FEWZ computation predicts the shape correctly, within scale uncertainties of order 
6-12%, where the default scale has the general feature of underestimating the relative abun¬ 
dance of high-ijx (>20 GeV) events at the 7% level. The shape in \y\ is well described by FEWZ. 

The absolute double differential cross section is presented in Table [2| The comparison with the 
NNLO computation of the FEWZ program is shown in Fig. [4] The scale uncertainties range from 
10-16% for qj > 20 GeV. The PDF uncertainties are of the order of 3% in the central rapidity 
region and decrease to about 1% in the forward region. The absolute cross section predicted by 
the NNLO program FEWZ agrees within the uncertainties with the measurement. 

A comparison of the measurements with the MadGraph and the POWHEG generators is shown 
in Fig. [5] The statistical uncertainties are smaller than the symbol size. The hatched bands rep¬ 
resent the systematic uncertainties of the measurement only. The two generators show opposite 
trends in qj. The MadGraph generator overestimates the data in the highest qj bins, whereas 
the POWHEG generator underestimates the data up to 20% in this region. Also shown are the 
absolute differential cross section predictions of MadGraph and POWHEG after normalising 
their inclusive cross-sections to the NNLO cross section by K factors, that are independent of 
qj and \y\. 


Table 1: Measured double differential fiducial cross section normalised to the inclusive fiducial 
cross section in units of GeV 
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8 Summary 

For Z bosons decaying to muons the double differential Z boson fiducial cross section in qj 
and | y\ has been measured in pp collisions at 8TeV. The results are compared to the next-to- 
next-to-leading-order predictions computed with the FEWZ program and they agree within the 
scale uncertainties. Deviations from the data of up to 20% at high transverse momentum are 
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Figure 3: The measured fiducial Z boson differential cross section, normalised to the inclusive 
fiducial cross section compared to the NNLO prediction of FEWZ. The first five plots show the 
cjj dependence in the five bins of \y\ and the last plot shows the tjj dependence integrated over 
\y\. The NNLO NNPDF23 PDF set with radiative corrections is used for the generation. We 
include data in qj up to 1 TeV, but have shortened the bin for presentation purposes. 
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Figure 4: The measured absolute fiducial Z boson differential cross section compared to the 
NNLO prediction of fewz. The first five plots show the i]j dependence in the five bins of \y\ 
and the last plot shows the qj dependence integrated over \y\. We include data in t]j up to 
1 TeV, but have shortened the bin for presentation purposes. 
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Figure 5: Normalised (left) and absolute (right) fiducial Z boson cross section, as a function 
of cjj, compared to predictions from MadGraph (red symbols) and POWHEG (blue symbols). 
MadGraph uses the LO CTEQ6L1 PDF set and POWHEG the NLO CT10 PDF set. The inclu¬ 
sive LO MadGraph and the inclusive NLO POWHEG cross sections are scaled to the inclusive 
NNLO cross section calculated with FEWZ by applying scale factors K^^ L q. 

observed in the MadGraph and POWHEG generators. The results are presented along with the 
full covariance matrix in order to enable their use in future fits of the PDF. The experimental 
uncertainties are significantly smaller than the current theoretical and PDF uncertainties. 
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